home *** CD-ROM | disk | FTP | other *** search
/ Info-Mac 4 / Info_Mac IV CD-ROM (Pacific HiTech Inc.)(August 1994).iso / Science / MAL ©P.f.Howden 1⁄1⁄89 / INTEGRATE < prev    next >
Encoding:
Text File  |  1989-07-12  |  1.0 KB  |  17 lines  |  [TEXT/MSBA]
  1. 1 CLS:PRINT"INTEGRATOR=Area (quadrature) or Line (arc) for single-value Functions.":INPUT "MENU: 0=A Function:1=Area under N DATA points:=";M:DEFDBL A-I,O-Z:IF M=0 THEN 4
  2. 2 CLEAR:DEFDBL A-I,O-Z:INPUT "ORDINATE SPACING H=";H:INPUT "N(ODD)=";N:DIM X(N):FOR J=1 TO N:PRINT "INPUT X";J;:INPUT "=";X(J):F=F+X(J)*(3+(-1)^J):NEXT J
  3. 3 PRINT "AREA INTEGRAL=";(F-X(1)-X(N))*H/3:PRINT"END :":PRINT:GOTO 2
  4. 4 PRINT "TYPE FUNCTION ON LINE 10.  eg: 10  F=X^3-2":INPUT "MENU: 0=AREA under curve: 1=LINE integral: ?";L
  5. 5 IF L=1 THEN INPUT "DIFFERENTIAL (.0001?)=";Z
  6. 6 GOTO 20
  7. 7 IF L=1 THEN X=X+Z:GOSUB 10:F1=F:X=X-2*Z:GOSUB 10:F2=F:X=X+Z:F=SQR(1+((F1-F2)/(2*Z))^2):RETURN
  8. 10 F=3*X^5-4*X^3
  9. 19 RETURN
  10. 20 INPUT "LOW LIMIT=";A:INPUT "TOP LIMIT=";B:INPUT "ACCURACY (.001?)=";D
  11. 30 H=(B-A)/2:X=A:GOSUB 7:P=F:X=B:GOSUB 7:P=H*(P+F):I1=3*P:N=1
  12. 40 S=0:FOR K=1 TO N:X=A+(2*K-1)*H:GOSUB 7:S=S+F:NEXT K
  13. 50 I0=P+4*H*S:IF ABS(I0-I1)<D THEN 70
  14. 60 P=(P+I0)/4:N=2*N:H=H/2:I1=I0:GOTO 40
  15. 70 PRINT "INTEGRAL=";I0/3:PRINT"END :":PRINT:GOTO 5
  16.                                        MACINTOSH LISTING
  17.